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Abstract 


The distortion of a two-dimensional bubble (or drop) in a corner of angle 6, due to 
the flow of an inviscid incompressible fluid around it, is examined theoretically^ The 
flow and the bubble shape are determined as functions of the angle the contact angle 

and the cavitation number Y^ The problem is formulated as an integrodifferential 
equation for the bubble surface* This equati^ generalises the integrodifferential 
equations derived by Vanden-Broeck and Keller^* The shape of the bubble is found 
approximately by using the slender body theory for bubbles presented by Vanden-Broeck and 

Keller^* When y reaches a critical value Tq( 6,6), opposite sides of the bubble touch 
each other* Two different families of solution for 7 < Tq are obtained* In the first 
family opposite sides touch at one point. In the second family contact is allowed along a 
segment* The methods used to calculate these two families are similar to the ones used by 

O A 

Vanden-Broeck and Keller and Vanden-Broeck • 

1. Introduction and formulation 


We consider the steady potential flow around a gas bubble or liquid drop in a corner of 
angle 6* The contact angle is denoted by B (see Figure 1). We shall write "bubble" to, 
mean either bubble or drop. We take into account the surface tension a at the interface, 
but we ignore the flow inside the bubble, assuming that the pressure is a constant 
throughout it* 
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Figure 1. Sketch of the bubble and the coordinates 


bubble is (x + iy)"'"# where a is a constant and x and y are Cartesian 


In order to formulate this problem we assume that the complex potential without the 
(x + ly) 

coordinates* 

6 

20 

We introduce dimensionless variables by choosing ( — as the unit length and 

»-6 pa^ 

20 

■jt[ — j) as the unit velocity. We also introduce the dimensionless potential bf and 
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8tr«am function Here, b > 0 ie a dimeneionXess constant to be chosen so that 

♦ ■ 4 and ♦ * • 4 stagnation points on the walls y » 0 and y ■ xtan6, 

respectively* We denote the streamline along the two ^Ali* and along the bubble boundary 

by ♦ * 0. In these variables b(f + !♦) ^ ~ (x + iy)^ at infinity or, equivalently 

7 1 

X ♦ iy ' (^) (♦ + i*)* (1) 

at infinity. 


The flow occupies the region # > 0 of the t plane, and the bubble boundary 
corresponds to the segment ^ 4 ^ ^ ^ axis ^ « 0. The problem of finding the 

flow consists of determining x + iy as an analytic function of ^ ^ in the half plane 

t > 0 satisfying Equation (1) at infinity. Then the bubble surface is given by setting 

♦ » 0 in x(t ♦ it) and y(t + it) and letting t range from ~ to -j. The contact 

angle conditions require that the bubble surface meets the walls at the angle 3, which 
yields 


y^ tan 3 as ♦ J 
*t -tan{e -6) as ♦ " 4 * 


( 2 ) 


On the bubble surface the pressure in the fluid, which is given by the Bernoulli 
equation, must differ from p^ by oX, where X is the curvature of the interface. This 
leads to the boundary condition 


P - ^ 
^s 2 


Pj, - ok on - ^ 


(3) 


Here, p^, P and q are, respectively, the stagnation pressure, the density and the speed 
of the fluid outside the bubble* In dimensionless variables (3) becomes 
2 


q“ « X - Y on " J ^ ^ 2# ♦ “ ^ » 

where Y is the cavitation number defined by 


« 

Pw P o 2*- 4 

*^b ^s r 2o>i 

o 

pa 


(4) 


(5) 


The problem can be further simplified by requiring the bubble to be symmetric about tho 

Y 

line y « X tan This implies that 

y^(4#0) ■ y^(-4,0)cosfi - x^(-^,0)sin5, o < ♦ < ^ , (6) 

By using Equation (6) we can restrict our analysis to the interval 0 < ^ < 4* 

2* Reformulation as an integrodif ferential equation 

It is convenient to reformulate the boundary value problem as an integrodif ferential 
equation by considering the function 

1-1 7 

(♦ + it) *(x^ + iy^) - {-5^) *1 , 

which is analytic in the half plane t > 0 and vanishes at infinity as a consequence of 
Equation (1)* Therefore, on t ” 0, its real part is the Hilbert transform of its 

imaginary part* The imaginary part vanishes on t * 0, Itl > 4 therefore the Hilbert 

transform yields 
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1-4 


♦ X^(^.O) - (-|) T ■ ? / «!♦' 

1-1 

1 0 (-♦') [-y.(^', 0 )co»« + ; .(f'« 0 )aln 6 ] 

+ i/ i i df . 

* -1/2 ♦ ♦ 

We now uea the eyminetry condition ( 6 ) to rewrite (7) in the form 

i « i 1 6 

.h •« 4 ■! e* ^ 1/2 1-4 , , 

X^(t»0) ■ (-y) -y ♦ + l_-_- / (^') y^( ♦' . 0) ( ’ ^ T + ^1 4, • 


( 7 ) 


( 8 ) 


Next 


we express the boundary condition (4) in terms of x^ and y^ noting that 
q2 ■ b^(x^ + Then (4) becomes 




*J + yj (xj + yj)^ 


I ♦! < Ts ♦ ■ 0 . 


(9) 


Now ( 8 ) and (9) together constitute a nonlinear integrodif fcrentiaX equation for Yx(^) 

1 ▼ 
in the interval 0< ♦«0, The contact angle conditions (2) complete the 

formulation of the problem for 

For t * »s the equation defined hY ( 8 ) (^) reduces to the integrodif ferential 

equation derived bY Vanden-BroecK and Keller^. The particular case ® ■ 5 represents half 
of a free bubble. 

For Y * S ■ 5 # the equations ( 8 ) and (9) reduce to the integrodif ferential equation 

* 2 

derived bY Vanden-BroecK and Keller'^. This case represents a quarter of a free bubble in a 
straining flow. 

The integrodif ferential equation defined bY ( 8 ) and (9) can be solved numericallY for 
arbitrarY values of y and 6 bY using the numerical procedures described bY Vanden- 
BroecK and Keller^' 

In the next sections# we shall find the shape of the bubble approximatelY by using the 
slender body theory for bubbles presented by Vanden-BroecK and Keller^. 

3. Slender body approximation 

For Y large the bubble tends to an arc of a circle of radius As y decreases 

numerical solutions show that the bubble elongates in the direction of the line which 
bisects the angle between the two walls. Then it develops a horn or spiKe which large 
curvature near its end* Finally when Y reaches a critical value Y^(0#Os opposite 

i V 

sides of the bubble touch each other. For ^ < y the contact point is at x«y» 0. For 
$ > j the contact point is away from x - y "0. Typical profiles for 6 ■ 1 and 
6 « 0 ■ y can be found in Vanden-BroecK and Keller^* These profiles were obtained by 
solving numerically the integrodif ferential equation of Section 2. 

For Y ^ Yq( 0*^) the bubble is slender. Therefore we shall use the slender body theory 

for bubbles presented by Vanden-BroecK and Keller^ to get an approximate description of the 

flow around the bubble. In the lowest order# the flow about a aynonetric slender bubble is 

approx irMi ted by the flow about a rigid plate lying along the center line of the bubble. In 

the present case the center line of the bubble consists of a straight segment of some 
length a lying along the line y " x tan6/2* We introduce the coordinates x'#y* (see 
Figure 1) and find the potential b4(x'#y') of the flow about these plates requiring that 

w 

at infinity b(^ t it) ^ (x ♦ iy)^. Evaluating the potential on the plate 
y* • 0, X* > 0 we obtain 
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( 10 ) 




ii 


2« 2f 1/2 

b4(x*,0) - j la”^ - x*"^) 


By dlfferantiating (10) wa find that tha flow spaad 
2 « 2 » 2 » - 1/2 

q(x’,0) - X* (a~^ - x*~^) . x’ > 0 . 


q on the plate is 


( 11 ) 


Before using q to get the bubble shape, we shall determine the length a. We do so by 
requiring the suction force F, exerted by the flow on the end of the spike, to balance 

2 

the surface tension 2o. As we see in Ref* 5 tp. 412, Eq (6.5.4)], F ■ ipA /4. Here 

1/2 6 

A, is the coefficient in the expansion b4 ^ Ar ' cos ^ in terms of polar coordinates 

with their origin at the end of the plate. Upon setting F » 2o and introducing 
dimensionless variables we obtain 

a 2 - i . (12) 

1 

From (10) we find A^ * — a , so (12) yields 

6 

j 2i-i 

a - (|) . (13) 

We next use (11) for q in (4) and approximate the curvature k by * x * ^ ^ * Here 

the equation of the bubble is y' * n(x*). Then (4) becomes 

4» 2t 

’'x'x* “ - x* - Tf . (14) 

At the end of the spike we require 

n(a) * 0 . (15) 


In addition the contact angle condition yields 
n* (a) « -tan(3 - -j) . 


(16) 


Here a is defined by the equation 

n(a) * otan ^ . (17) 

The function n(x') is easily obtained by integrating (14) twice with the auxiliary 
conditions (15) and (16). In the particular case 6 » w, the result of the integration is 

n(x‘) * (a^ - x*^)(y - l)/2 ““a(a + x*)log(a + x'i 

- i a(a - x')log(a - x‘) + a^log2a + (a - x*)tan(0 - j) , (18) 


For 6 ■ w, (13) becomes 


a » — . (19) 

Vanden-Broeck and Keller^ have shown that the approximate solution (17), (18) is in fair 
agreement with the exact numerical solution of (8) and (9) for ^ ^ “J * 7 ^ 7 q(6,»). 

For Y < Yq( 3,5) (14) -(16) yield unphysical profiles in which opposite sides of the 

bubble cross over. In the next two sections we construct physically acceptable families of 
solutions for Y < Yq(3,3). We shall present these results in the important particular 

case 6 > w. 


4. Solution with one point of contact 

To obtain solutions for y < Yq(*, 3) we require the free surface to be in contact with 

itself at one point. Then the bubble contains a small sub->bubble near its tip (see Figure 

2). We denote by c the x* coordinate of the contact point. 

We describe the profile of the bubble by the equations y* ■ nj^(x*) 0 < x* < c and 

y’^n^ix*) c<x*<a. Then by symmetry we have 
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Figure 2. Profile of the bubble with one point of contact for y • Yq ^ -1.7 and 
6 » 2«/3. The vertical scale is the same as the horizontar ccale. The 
cavitation number in the sub-bubble is equal to Yq* 



Figure 3. Profile of the bubble with one point of contact for Y » -4.5 

and 8 ■ 2w/3. The vertical scale is the same as the horicontal 
scale. The cavitation number in the sub-bubble is u ■ -0.6. 
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( 20 ) 

( 21 ) 


ni(c) - r\^(c) « 0 

n|(c) - 12 ( 0 ) >■ 0 . 

The conditions (15) and (16) yield 

H2(a) - 0 (22) 

n*(0) - -tan(B - |) . (23) 

The functions n^(x’) and n 2 (x') are obtained by integrating (14) twice* The four 

constant oL integrations and the value of c have to be evaluated to satisfy the six 
conditions (20^-(23)* This yields a system of six equations with five unknowns. Therefore 
we cannot expect this system of equations to have a solution for any value of Y other 
than Yq{ 0 ,«). 

The physical reason why the problem does not have a solution for Y ^ Yq is that it 
requirea the cavitation number in the sub-bubble to be the same as in the main bubble. It 
is to be expected that the cavitation number within the sub-bubble will have some value 
other than y» which we cannot prescribe. Following the general philosophy of the method 

used by Vanden-Broeck and Keller^ we shall introduce the unknown cavitation number u in 
the interval c < x* < a. 

Integrating (14) twice we obtain 

n^(x*) - (a^ - x'^) 1 a(a + x')log(a x‘) 

- i a(a - x')log(a - x‘) + A + Bx' , (24) 

H 2 (x') « (a^ - x*^) 2 ^ i a(a + x')log(a + x‘) 

- j a(a - x')log(a - x') + E + Dx* . (25) 

Here A, B, E and 0 are the four constants of integration. Using the six conditions 

(20)-(23) we obtain a system of six algebraic equations for the six unknowns A, B, E, 

D, w and c. This system can easily be solved and yields a unique solution for any y in 

2 X 

the interval -• < ^ < Yq( 0,«). Typical profiles for ® “ 3 ^ shown in Figures 2 and 

3. The value of Yq is approximately equal to -1.7. As y decreases the sixe of the 

sub-bubble increases and the size of the main bubble decreases. For y ■ -**« M « -0.39 
and the main bubble vanishes. It is interesting to note that the present solution also 
exists in the interval ^ ^ ^ Here Y* is the value of y for which 

A similar result was found by Vanden-Broeck and Keller^. 

The results are summarized in Figure 4. The solution before contact described in 

Sections 2 and 3 correspond to the interval < Y < •• It is represented by the straight 

line M ■ Y in Figure 4. The othtjr curve in Figure 4 corresponds to the present 

solution. It exists in the interval < y < Y*. Therefore there are two possible 
solutions ir. the interval < Y < Y*. 

5. Solution with an interval of contact 

In this section we derive another solution for Y < Yq( 6 ,x) by requiring the bubble to 
be collapsed between x* * f and x* ■ g (see Figure 5); we describe the profile of the 
bubble by the equations y* ■ nj^(x'), 0 < x* < f and y* * 1 ^ 2 (x*), g < x* < a. The 
functions nj^(x') and n 2 (x*) '^^must satisfy the following conditions 


H2(a) ■ 0 (26) 
n|(0) - -tan(6 - j) , (27) 
n^(f) - m 0 , (28) 
n^(f) - n^(q) - 0 . (29) 
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Fiqure 4, The cavitation number u as a function of y. 



Figure 5. Profile of the bubble with one segment of contact for y ■ *3.0 and 
0 * 2«/3. The vertical scale is the same as the horisontal scale* 

The values of f and g are respectively 0*19 and 0.47* 

The functions n^(x*) and obtained by integrating (14) twice. They are 

therefore given by the relations (24) and (25). The six constants A, B, E, D, g and f 
are found by satisfying the six conditions (26)-(29). We note that the present solution 
can be found with the same cavitation number everywhere. 

2 s 

A typical profile for ® shown in Figure 5. As y decreases the sizes of the 

main bubble and of the sub-bubble decrease. Furthermore the length of the contact segment 
increases as y decreases. For y * the bubble reduces to a straight segment of 

length a lying on the x* axis. 
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Finally l«t ua mention that tha equilibrium of forces require the segment of contact to 
be a **film of impurities** characterised oy a surface tension equal to 2o. This is very 
unlikely to occur in reality* Therefore the bubble with a segment of contact is physically 
unrealistic* However# this mathematical solution is physically relevant to describe the 

deformation of an inflated membrane* For details see Vanden-BroecK^* 
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